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The electrical conductivity is calculated for regular inhomogeneous two component isotropic 
medium in which droplets of one phase with conductivity <T2 are embedded in another, with con- 
ductivity 01. An expression is formulated which can be used in many different situations, and is 
of particular relevance in the case where the relative proportion of the components is temperature 
dependent and varies over a wide range. The behavior of effective conductivity depends on the 
spatial arrangements and the shape of the inclusions. 



I. INTRODUCTION 

The determination of the effective conductivity a e ff of spatially inhomogeneous heterophase systems is an old, yet 
increasingly important problem of theoretical physics. With the advent of new nanoscale probes of condensed matter 
systems, it has been apparent that many very diverse systems which were previously thought to be homogeneous are 
in fact either statically or dynamically inhomogeneous. The effective conductivity in such cases cannot be dealt with 
in terms of homogeneous medium theory, is not trivial, and solutions are known so far only in some rather special 
cases. Different aspects of the theory and different limiting cases are extensively discussed in the Rcf. [1]. 

In this paper we focus on the problem of calculating the effective conductivity of an inhomogeneous two-dimensional 
(2D) plane. The classical problem could be formulated as follows. Let us assume a 2D system which contains a mixture 
of N (N > 2) different phases or materials with different conductivities <7j, i — 1,2, ...,iV. The arrangements of 
different phases could be random or regular. The question which we wish to address is how the effective conductivity 
of the plane depends on the conductivities of the phases, their concentration and the spatial arrangements of the 
phases. 

In the past, a number of different approaches have been used to tackle this problem. The exact result for the effective 
conductivity of a two component system with a symmetric and isotropic distribution of components was obtained by 
Dykhne [2] . He found that effective conductivity of the system is determined by the simple relation o~ e f / = \/(J\<Ji- A 
"symmetric" distribution, applied to this problem is one in which the two components can be interchanged without 
changing the end result. Obviously, one requirement for a symmetric distribution is that the two components have 
equal proportions, but it also means that more general cases cannot be considered with this model. 

Further investigations have shown that more general duality relation is valid for 2D heterogeneous conductors 
then that initially considered by Keller and Dykhne [1]. More recently it was shown that more general relation for 
the tensor of the effective conductivity exists which is valid for multicomponcnt and anisotropic systems [3,4]. The 
effective conductivity of several examples of ordered two-component systems was also calculated exactly [5-7] . It was 
shown [5,6] that for a chess-board plane and for a plane constructed of triangles, the relation derived by Dykhne is 
also valid. 

A similar relation to the Dykhne formula for the effective conductivity of a system consisting of randomly distributed 
metallic and dielectric regions near a metal-to-insulator transition was derived by Efros and Shklovskii [8]. They 
generalized the expression of Dykhne on the basis of scaling arguments to the case of arbitrary concentrations of the 
two phases near a percolation threshold, so that the effective conductivity becomes: 

0eff = cti(ct 2 /cti) s , (1) 

where s is a universal scaling exponent. Critical exponents are relatively well known as well for this type of system 
[9]. This relation is not applicable when the system is driven away from the percolation threshold and the general 
solution of the effective conductivity of an inhomogeneous medium thus remains an open problem. 

Sen and Torquato [10] derived an expression which allows for explicit calculation of the effective conductivity 
tensor from n-point probability functions 5„(r l7 . . . , r„). These functions give the probability of points at ri, . . . , r„ 
belonging to the same phase, and are thus uniquely determined by spatial distribution of the phases. Unfortunately, 
the application of this method is limited as the computations involving n > 5 are fairly time consuming. 
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Different expansions of the effective conductivity in terms of small parameter have been used in the past [1,11,12]. 
In most of the cases the low order terms weakly depend on microgeometry. A diagrammatic expansion for the effective 
conductivity was developed by Khalatnikov and Kamenshchik [13] which promises to give more generally applicable 
results. The perturbative approach seems to be quite effective since it allows us to analyze random and nonsymmetric 
distributions with different conductivity. 

The problem was discussed for the case N = 2 and N — 3 on the basis of numerical calculations [3,14] as well. It 
was shown that the effective conductivity for N — 3 is not universal and depends on the spatial arrangements of the 
phases. We have employed boundary element method for efficient numerical treatment of two-dimensional multi-phase 
systems with arbitrary arrangement of phases. More details on the method and its results can be found in [15]. 

In this paper we consider the conductivity of a two-phase system in two dimensions for a wide range of concentration 
and conductivities. One phase is assumed to be composed of droplets (of different shapes) with conductivity 02 
embedded within a medium of conductivity u\ (see Fig. 1). We begin by calculating the effective conductivity e // 
using a perturbation theory approach with the two phases having volume fractions (1 — v) and v respectively. Since 
the problem is linear, we can introduce a dimensionless conductivity 0, measured in units of o\ = 1, and the effective 
conductivity a e ff is a function of a = 02/01 = 02 and v. The volume-averaged conductivity = y J crdV is given 
by: 

a = (l-v) + va (2) 

Assuming that the conductivity of the two phases is not vastly different \a — 1| <C 1, the effective conductivity 
could be calculated by perturbation theory [13]. To apply perturbation theory we rewrite the spatial dependence of 
conductivity as: 

a(r) = a(l - a(r)), (3) 
where a(r) = £Mz£. Then, assuming that the spatial distribution of conductivity is uncorrelated, we obtain: 



/ 



dva(r)a(v + r') = ^ v) 6(r') (4) 



After a straightforward calculation, the conductivity up to second order in a is given by: 

(q-i)Mi-^) m 

This result has been known for many years and was derived for the dielectric function of dielectric mixtures [16]. In 
Refs. [11,12] it was also derived using a systematic perturbation expansion which showed it to be exact to second 
order in a. The second term in Eq.(5) represents the first non- vanishing contribution due to the inhomogeneity of 
the distribution of the phases. For the case v = 0.5, the result coinsides with the expansion of the exact expression 
for the conductivity up to the second order in (cr 2 — 01) [2]: 

Veff = (6) 



II. CONDUCTIVITY OF A REGULAR ISOTROPIC TWO-COMPONENT SYSTEM IN 2D. 

Next, we calculate exactly the effective conductivity of the plane with different regular isotropic distributions. As 
before, let us consider a 2D plane which is constructed from two different phases with different conductivities 01 = 1 
and 02 = 0- The regions with conductivity 02 have a circular shape with radius R and form a regular square lattice 
with the period a as shown in Fig. la. Changing the radius R from to R = a/2 we can change the volume fraction 
of the second phase from v = to a critical concentration v c = 0.785 whereafter the regions with conductivity 02 
start to overlap and a percolation threshold is reached. In the case of metallic droplets, the total charge density 
should be zero, while a finite charge density can accumulate on the surface between different phases. This allows us 
to formulate the integral equation for the surface charge density [5,6]. Let us define a surface charge density by the 
relation p{9)Rd9 = dp(9), where dp(9) is the charge on the small part of the surface between the two components 
with the length dl = Rd9. Taking into account that the scalar potential at the point r is determined by the relation 
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cj) = E x-2 j d 2 r ln|r-r'|p(r') (7) 

where In |r — r \/2-k is the 2D Green's function. The boundary conditions on the surface between two phases are [16]: 

El - El = 4np(9) (8) 

a x El = a 2 El (9) 

Substituting r = (ma + R cos (0 ))i+ (na + Rsin(9 ))j and r = i?cos (9)i + Rsin to Eqs.(7-9) we obtain integral 
equation for the surface charge density in the following form: 



p(9) = — [E cos (9) + 2r ]T / d9 p(9 ) 

n in = — rv-i 1" 



■.m + r(cos (9') — cos (9)) + i(n + r(sin (9') — sin (9))), 



(10) 



where r = R/a, and k = . As it is shown in the Appendix the sum over m can be calculated exactly and the 
integral equation for surface charge density will be reduced to the following form: 

p(9) = — [E cos (9) + 2r V / d6 K(n,9,0 )p{9 )] (11) 

where 

/ _ cos(0)sin(27rr(cos(0') - cos(0))) + sin (0) sinh (27r(n + r(sin (6>') - sin(0)))) 

(n ' ' cosh (27r(n + r(sin [9') - sin (9)))) - cos (27rr(cos (0') - cos (9))) (1 ' 

Expanding the surface density p(9) in terms of Legendre polynomials Pi (cos (9)), and taking into account that 
p(-9) = p(9) and p(ir -9) = -p(9): 

OO 

pW=^c 2; _ 1 P 2i _ 1 (cosW) (13) 
i=i 

we obtain the following linear set of algebraic equations for the coefficients c%-\ 

|^ = A [ 2^ 1+2r g C2fe _ 1 ^ fe] (14) 

fe=i 

where 

K ttk = V / d9 / dflA-(n,e,e')sin^ I _i(cos(e))i^ fc _i(cos(fl')) (15) 

Solving Eq.(14) taking into account a finite number of Legendre polynomials we obtain surface charge density 
Eq.(13). As a result the effective conductivity is evaluated by calculating the total current j = u\E n = E n through 
the semicircular surface with the radius R — a/2 (see Fig. la). Calculations, similar to that of Eq.(12), lead to the 
expression for effective conductivity: 

a eff = — / d9[cos (9) + — / d9K («. 6 W )] ( 16 ) 

/ / cos(6»)sin(2 7 r(rcos(6'')-cos(0)/2)) + sin(6')sinh(27r(n + rsin(6'')-sin(6')/2)) 

(n ' ' ' ~ ^ cosh (2tt(u + r sin (9' ) - sin (9) /2)) - cos (2n(r cos (9') - cos (9) /2)) ( ' 

The result above applies to the case of a uniform distribution of circular droplets within the plane. To see how the 
effective conductivity depends on the shape of the regions with conductivity <T2, we have performed calculations for 
the case where of circular droplets we substituted with squares, triangles and rhombuses with the ratio of diagonals 
tana = a/b where a and b translation vectors along x and y respectively (see Fig.lb,c,d). In all these cases Eqs.(10- 
17) are slightly modified since in polar coordinate system r(9) is a function of angle. Contrary to the case of circles, 
percolation threshold for cases b,c, and d is v c = 0.5. Note that in the case of rhombuses the lattice is anisotropic 
and al} f ^af ff . 



where 
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III. DISCUSSION. 



The results of the calculations of the effective conductivity are presented in Fig. 2 as a function of a v for different 
values of volume fraction v. It is easy to check that the results satisfy the generalized duality relation [3,4]: 

^((71,(72)^(1/(71,1/^) = I- (18) 

For the cases of circles, squares and triangles a 22 — a 11 . In the case of rhombuses a 22 {a) = a 11 (ir/2 — a). Fig. 2 
(a,b,c,d) shows that for small k, perturbation theory [11-13] (Eq.5) gives the correct result independent of geometry. 



A. Approximate expression for the effective conductivity. 

Although the predictions in Fig. 2 represent the results of a precise numerical calculation, they are not very tractable 
when it comes to comparing with experimental data, being the result of numerical calculations. It is therefore helpful 
to try and obtain a functional form for describing the behavior predicted in Figure 2, which also includes all the 
relevant parameters, such as volume fraction v, and the two conductivities <j\ and oi- Such an expression can then 
be used for a wide range of problems, provided the range of validity is taken into account. We describe the properties 
of such a heuristically determined function and determine its range of validity in terms of parameters v, o\ and 02 ■ 

As can be seen from the Figure 2, the dependence of the effective conductivity on a shows similar behavior, 
independent on the particular geometry of the phases. First, we observe that when k is small all the curves are 
linear in o v with the same slope. In the relatively wide interval of sigma (0.1 < a < 10) the effective conductivity is 
determined by the equation: 

<Teff(*) = <T[ 1 - V) a$. (19) 

The range of applicability of this formula becomes wider as we approach the percolation threshold v c . When a = 
02 /01 3> 1 effective conductivity saturates at a sat ■ cr sa t is not universal and depends on the geometry. Recently it was 
pointed out that for the case of circles for v < 0.5 in the whole range of a effective conductivity may be approximated 
by the formula [18]: 

CT e //(«) = (1 -vk)/(1 + vk). (20) 

To derive approximate expression for effective conductivity we assume that Eq.(20) remains correct if we substitute 
instead of v the effective volume fraction v e ff(n, v). We require that v e ff{n,v) w v for k — > or v — > 0, and 
v e ff[K, v) rts ^ \+°J c f° r v ~ > v c to satisfy Eq.(19) which is valid at v = v c . It is easy to see that the function: 

1 _ f i-ti-pMia 1 '' 

Veff{K,v) = v+- — — — „ e - {l-p{v))v c , (21) 

1 + {l+(l-p(v))K ) 

where p{v) — > as v — > v c and p{v) — > 1 as v — > 0, satisfies all of the above requirements. Function p(y) is not 
universal and depends on the geometric shape of the region with conductivity a and on the particular arrangement 
of these inclusions in the 2D plane. In Fig. 3 we plot p{v) as function of 1 — vjv c for the cases a,b,and c respectively. 
The case d is different, since the effective conductivity is anisotropic. As it is clearly seen from Figure 3 the behavior 
of the function p{v) for circles (case a) is different from the cases of squares and triangles (b,c). On the other hand 
in the cases b and c p(y) shows similar behavior. 



B. Shape dependence of the effective conductivity. 

The function p(y) is connected with the value of a sat — (1 + v e ff(K = — v e ff(K = — 1, v)). Therefore the 

behavior of the function p{v) close to percolation threshold should be different for different geometries. In Fig. 4 we 
plot the value of a sat as a function of (1 — v/v c ) for the case of circles, squares and triangles. There is an important 
difference between these two cases. In the case of circles a sat diverges as a power of (1 — v/i> c )~ k (k w 0.5). For the case 
of squares and triangles this behavior is logarithmic. In both cases close to percolation threshold a sat is proportional 
to average inverse distance between boundaries of the neighboring circles or squares <j sat cx J dy/(l — 2f(y)) and 



4 



f(y) — yj (r 2 — y 2 ) for circles and f(y) = r — \y\. Here we assume that the period of the system is 1, and the 
dimcnsionlcss size of the circle and the square is r. Direct integration leads to the following results: 

a sat oc (tt/2 - arccos((l - v/v c ) l l 2 )/{l - v/vc) 1 ' 2 - tt/2 (22) 

for the case of circles and 

a sat cx-ln(l -(!//!/<=) V 2 ) (23) 

for the case of squares (Fig. 4). For the case of triangles (c) the asymptotic is similar to Eq.(23) with different numeric 
coefficients. Interestingly this observation suggests that behavior of the function p{v) is different depending on the 
curvature of the embedded regions. 



IV. CONCLUSION 



From calculations of the effective conductivity of inhomogeneous two-phase systems in two dimensions we find 
that the results of precise numerical calculations can be approximated by a universal function for <r e ff Eqs. (20,21), 
where the function p{v) depends on the spatial arrangements of the 2D plane and on the shape of the inclusions with 
conductivity a. It is shown that in a large interval of conductivity a, the effective conductivity cr e // is determined by 
the spatial average of the logarithm of individual conductivities. The closer the system is to the percolation threshold, 
the larger the range of validity of this result. For large values of conductivity cr, cr e // saturates at the value a sat - 
The value of a sat near percolation threshold is determined by the average inverse distance between boundaries of 
neighboring regions with conductivity a in the direction of the field (Eqs. 22, 23) 

The model which we have developed is quite generally applicable and can be applied in some interesting situations, 
such as cuprates and other two-dimensional complex transition metal oxides which exist near a phase-separation 
threshold. Importantly, there appears to be a significant amount of experimental evidence that many anomalous 
properties of oxides are associated with the coexistence of two or more phases. The application of the presented 
model may help in understanding the transport properties of such systems. 



V. APPENDIX 



Here we show how the sum over m in the Eq.(10) could be calculate exactly. Let us represent the sum in the form: 

^ imn. A- ft A- irv\ ' 



m— — oo 



where (3 — r(cos(6* ) — cos(#)), a = n + r(s'm(9 ) — sin(6>)). Sum over m is calculated using the definition of di-gamma 
function. As a result we obtain expression for the sum: 

S = $t(exv{i6)[ip(-/3 - ia) -ip(l + /3 + ia)}) = ir$i(exp(i6) cot(ir{/3 + ia))) 

Calculating imaginary part of previous equation we arrive to the results Eq.(12): 

_ ^cos (9) sin (2-jt/3) + sin (9) sinh (2ira) 
cosh (2na) — cos (2n(3) 
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Figure 1. Spatial arrangements of phases with conductivities ai = 1 and oi = a for four considered cases. 
Figure 2. Effective conductivity of the plane as a function of <j v for different volume fractions and four considered 
geometries. 

Figure 3. Dependence of the function p(v) on 1 — v /v c for the cases a,b and c. 

Figure 4. Saturated effective conductivity as <r — > oo for the case of a, b and c. Full, dotted and dashed lines show 
different analytical asymptotical behavior for these cases. 
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